Magnon mechanism of Josephson coupling in SFS structures 
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It is shown that Josephson coupling in SFS junction due to electron-magnon interaction remains 
at a distance, when the usual proximity effect decreases exponentially. We obtain expression for the 
Josephson energy, which contain the parameters of the magnon spectrum and allow to estimate the 
value of the maximum superconducting current. 
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The expected physics of the proximity effect in the 
structures of S-superconductor-ferromagnetic metal pj is 
based on the fact that the wave function of a Cooper pair 
in a metal oscillate with the distance from the boundary 
with the superconductor. This is the analogue of the be- 
havior in the superconducting LOFF state 0, 0] ■ The 
presence of impurity scattering in the metal leads to an 
exponential decrease of the wave function of Cooper pair 
on the oscillation length, equal to = y/D/h, where D 
is the diffusion coefficient, and h is the energy of the fer- 
romagnetic splitting. Experimental investigation of this 
pattern is well established see also the review j^-Q • 

Recent experiments on SF structures found [8l4l0l] , that 
the Cooper pairs in a superconductor-magnetic metal 
penetrate to a distance much greater than the length 
Lh- The explanation is that in these structures magnetic 
state is characterized by a noncollinear ordering. In Q 
it is a helical ordering, in [l(| it is artificially created in 
magnetic multilayers. In this case, odd frequency triplet 
state penertates into magnetic metal 0, [ll[ . The spatial 
symmetry of this state determines its insensitivity to a 
potential scattering. 

Mesoscopic fluctuations of Josephson current are pre- 
served on all scales associated with both the elastic in- 
teractions and ferromagnetic splitting. Therefore, in SFS 
junction, in which distribution of the of the supercon- 
ducting phase difference is correlated with the mesoscopic 
fluctuations, there is an average Josephson current [12| 
even at thicknesses d much larger than Lh- 

In this paper we study the influence of electron- 
magnon interaction on the Josephson energy of the SFS 
junction with thickness d >> Lh- Usually, when con- 
sidering the problem of the SFS contact, ferromagnet is 
modeled as a metal with a built-in spin splitting. This 



approximation neglects the fact that due to electron- 
magnon interaction electronic states with a given spin 
projection are not eigenstates. Therefore, the wave func- 
tion of Cooper pairs penetrating into the ferromagnet 
from S-superconductor contains a component, which does 
not oscillate in the exchange field and, as it is shown 
by calculation, does not decreases exponentially over a 
length Lh- Here we calculate corrections to the ther- 
modynamic potential of SFS junction due to electron- 
magnon interaction. 

We consider the s-d model for a ferromagnet with 
Hamiltonian 

H sd = J J dr^+(r)a lJ V^(r)S(r) (1) 

Here a are the Pauli matrices, S(r) = J^k — r i)^k is 
the density operator of spins, localized at the points r^. 

We assume ferromagnetic ordering of localized spins in 
the direction z and corresponding splitting of the electron 
spectrum of ferromagnetic h = J(S). (S) is the density 
of localized spins. 

The transverse part of the Hamiltonian ([T]) with the 
help of Nambu operators 



might be written as 




(3) 

The second-order perturbation correction to the thermo- 
dynamic potential has the form 
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The integration is over the ferromagnetic contact re- Here -Dj,fc(Ti2,ri, r 2 ) is the Matsubara Green's function 
gion |x|< d/2. We assume that it is in the y, z plane. 
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FIG. 1. The diagram corresponding to magnon contribution 
to the thermodynamic potential, averaged over the random 
scattering potential. The wavy line in the figure denotes 
the magnon propagator. The solid lines correspond to 2 x 2 
electron Green's function, averaged over the random poten- 
tial. Dotted lines correspond to the impurity scattering. It is 
shown element of the diffusion propagator. 



for the spin operators 

A,fc(na,ri > r 2 ) = -{T T S*(ri, ri)5 fe (r 2 , r 2 )} (5) 
Electron Green's function, defined as 

Gi(Ti2,ri,r 2 ) = -{IVtfitri.riJtf+fa.ra)} (6) 

is the 2x2 matrix. 

Operators ^i(r) and \& 2 (r) differ only in the direc- 
tion of the electron spins, therefore G\ (r, ri , r 2 , h) — 
G 2 (t, n,T2,-h) 

The figure shows a diagram corresponding to the ther- 
modynamic potential averaged over the random po- 
tential. We apply the standard technique of averaging 
described in the review 0^3 and the references cited 
therein. 

The averaging is performed over the scattering poten- 
tial, which in the Nambu representation is <j z V(y). It is 
assumed that (V(r)V(r')) = 27r ^ QT ^(r — r'), where r and 
Vq are the mean free path and density of states per spin 
at the Fermi level. 

The equation for the electron Green's function, aver- 
aged over the scattering potential has the form 

(ie n -h-H a z + ^g 1 (r))(G 1 (r,r',e n ))=S(r-r') (7) 

Here, e n = (2n + 1)ttT is the Matsubara frequency, T 
is the temperature, and Hq is the Hamiltonian of free 
electrons. 

The off-diagonal elements of the Green's function de- 
scribe the penetration of Cooper pairs in the normal 
metal over the length of the Lh = y/D/h. In our case, the 
thickness of the ferromagnetic of metal d >> Lh, ther- 
fore both superconductor-ferromagnet interfaces might 
be a considered independently. 

Near the boundary with the superconductor,which or- 
der parameter has a phase ip, the matrix g\ can be written 



9l (r) = Gi r,r,e„ = ./'^ J \ 



(8) 



The matrix elements are related by g\ + f± = 1. 

At a distance from the boundary larger than the Lh, or 
in the case of a weak proximity effect, the matrix elements 
are |/i|<< 1, gi = sign(e n ). 

The function /i(r,e n ) at |e„"r|<< I hr < 1 might be 
determined from the Usadel equation, which in the case 
of a weak proximity effect has form 

{-D\/ 2 + 2(e„ + ih)sign{e n ))h{v, e n ) = (9) 

Here D is the electron diffusion coefficient. Note that 
/2(r, e n ) might be obtained from /i(r, e n ) by changing 
sign of h —> —h. 

In the case of low transparency of the SF interface, 
the boundary conditions at the contact between the 
superconductor-ferromagnetic metal is fl3j 



7nV/i(r,e„) 



V^ 2 



(10) 



Here (7)" 1 is the ratio of the resistivity of the ferromag- 
net to the resistance per unit area of surface, n is the 
normal to the surface. A is the modulus of the supercon- 
ducting order parameter. 

In a bulk of ferromagnetic metal the diffusion propa- 
gators appearing in susceptibility (see figure) 

11= (Sp(CT z Gi(ri,r 2 ,e„)er 2 G 2 (r 2 ,ri, e n >))), (11) 

which are proportional to (1 ± cri 1 '' ) (1 =p tri 2 ') do not con- 
tain h. Accordingly, only those contributions are not 
damped at the length Lh and should be considered for 
the calculation of 5f2 for junction with d » Lh- In this 
case the frequencies e n and e n / must satisfy the condition 
e„e n / > 0. 

<7i(r) is slow varying function of coordinates at the 
mean free path distance. Neglecting it's gradients we 
obtain for vertex part of (|11[) expression 



d 3 p 



;(Gi(r,p,e„))cr z (G 2 (r,p,e n /)) = ttv t<t z (l-g^) 

(12) 

Substituting in (fl2|) definition ((5J in the case of a weak 
proximity effect we obtain for e n t n ' > 

/ e iv \ 

o- z (l-5l(r)32(r)) ~ -i \ ilp Q j (/i(r, e n )-/ 2 (r, e„0) 

(13) 

Thus, the entering into long-range contribution to the 
II vertex parts are different from zero only near the sur- 
faces of contact. 

The equation for the diffusion propagator, neglecting 
Andreev reflection (51(f) = sign(e n )a z ) is obtained using 
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the Green function (J7J with e n e n / > 0. It has the form 



47T1/0 7 



(-£)V 2 +|n n , n ,|)P(r,r / ) O„, n = ^(r-r') (14) 

Here, Qn,n' — £n' • 



In the limit d >> Lh = y/D/h in the integration over 
the spatial coordinates in Q slowly varying function of 
the coordinates D and P can be set equal to their values 
in the x = id/2. The integration of vertex parts over x 
is reduced to the replacement 

(/i(ri,e n ) - / 2 (ri,e n /))(/ 2 (r 2 ,e n <) - /i(r 2 ,e„)) -> 

\ \ 2 



(15) 




There are two contributions to the thermodynamic po- 
tential. 

If r"i and r 2 are on opposite surfaces of junction, than 
the phase dependent susceptibility has the form 



<5iil(ri,r 2 



D 



2hj \y/^+ 
xP(ri,r 2 ,n„ in ')cosiy9i 2 



(16) 



When the coordinates ri and r 2 belong to the same 
surface, phase dependent contribution arises after taking 
into account Andreev reflection from the opposite sur- 
face. Graphically, this means inserting a Hikami block 
containing anomalous part of the Green's functions in 
the diffusion propagator shown in Fig. Rate of the re- 
flection is DV/i(r)V/ 2 (r). As a result, 



(5 2 i7(ri,r 2 ) 




cos2(^ 12 J drP(r 1 ,r,f2 Bin /)V/ 1 (r)V/ a (r)P(r, r 2 ,ft„, n ,) 

(17) 



The integration over r is near surface, which is opposite 
to that of the ri and r 2 . Here again we can put the diffu- 
sion propagators are equal to their values at the surface 
and integrate DV/i(r)V/ 2 (r) over x. 

Since when h >>|e„|, |e n /| <5i 77 and (S 2 77 are even func- 
tions h, the thermodynamic potential depends on a com- 
bination of spin Green's functions as 



Dx,x(q, e« - £«' ) +D v<y (q, e n - e n > ) 



2(S)E(q) 



E{q) 2 + (e r 



(18) 

Here E(q) = E a + D s q 2 is the magnon energy, E a takes 
into account the anisotropy energy and the external mag- 
netic field, (S) is the spin density. 

Let consider the limits of 1). weak anisotropy 
\J D s /E a > Lh, when magnon propagator is slowly vary- 
ing on the length Lh, and 2). temperature, such that 
T>D s /d 2 ,E a . 

In this case summation over frequencies in thermody- 
namic potential Q might be restricted by terms with 

Normally D >> D s , so in this limit there might be 
any possible relationship between the thickness d and the 
coherence length y/D/T. 

After substituting in (|4]) expressions (fT6"| , (fl7|) and the 
magnon propagator, calculated for an infinite medium, 
and integrating over the SF surfaces, we finally obtain 



5Q = 5^ + 52^2, 



(19) 



where per unit area 



8if2 _ Dv cos (fi2 ^ A 2 



S 8j 2 D s (S) ^A 2 + e 2 n 




(20) 



and 



s 2 n 
s 



Dv cos 2<pi2 I DD 
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2e n \ 
D 



d\ (21) 



In deriving these expressions we used the relation h = 
J(S) between the sd interaction constant and the value 
of ferromagnetic splitting. Integration of diffusion and 
magnon propagators over surfaces is reduced to integra- 
tion over t. 

Note that minimum of corresponds to tp = tt state, 
and minimum is achieved at tp = 0, 7r. 

The spin-orbit interaction has two consequences. The 
presence of a gap in the magnon spectrum, corresponding 
to the anisotropy energy E a for d > y/D s /E a leads to 
the factor exp(— dy/ E a /D s ) in the expression (|20)) . The 
spin-orbit scattering of conduction electrons can be ac- 
counted for by the substitution |2e„|— >|2e„|+T so 1 in the 
exponents in (|20p and (|21[) . t so is the spin relaxation 
time of conduction electrons due to spin-orbit scattering. 



When d is greater than the length of the spin relaxation 
of conduction electrons L so — \JDt so the first contribu- 
tion decreases as exp(— d/L so ). In this case, the second 
contribution decreases with increasing spin-orbit scatter- 
ing faster than the first as exp(— 2d/L so ). 

Let estimate the value of the maximum superconduct- 
ing current corresponding to the expressions (|20[) and 
(12"T1) with y/D/T, L so > d, A > T, when the sum over 
the frequencies gives a contribution of order of unity. 

At specific resistance of the ferromagnet ~ 0.1/xxcm, 
the factor Dvq is ~ 2 x 10 10 cm _1 . For 3d metals, magnon 
spectra have D s ~ 10~ 17 eV cm 2 |l4j. At spin den- 
sity (S) ~ 10 22 cm" 3 and temperature ~ 1°K we have 
8d"(S) ~ 2°^- 7 is the ratio of the mean free path 
in the ferromagnet to the boundary transmission coeffi- 
cient [l3|. If 7 ~ 10~ 5 cm than for area S ~ 10" 8 cm 2 
maximum superconducting current is a few /iA. 

Note that because of the smallness of the D s the length 
yjD s /E a may be small even for weak anisotropy, or mag- 
netic field. In this case, contribution to (I2U1) is less than 
the contribution of (|2Tj) . 
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